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Abstract: We compare, using a three neutrino analysis, the allowed neutrino os- 
cillation parameters and solar neutrino fluxes determined by the experimental data 
available Before and After Neutrino 2004. New data available after Neutrino 2004 in- 
clude refined KamLAND and gallium measurements. We use six different approaches 
to analyzing the KamLAND data. We present detailed results using all the avail- 
able neutrino and anti- neutrino data for Am 2 ^, tan 2 # 12 , sin 2 #13, and sin 2 rj (sterile 
fraction). Using the same complete data sets, we also present Before and After de- 
terminations of all the solar neutrino fluxes (which are treated as free parameters), 
an upper limit to the luminosity fraction associated with CNO neutrinos, and the 
predicted rate for a 7 Be solar neutrino experiment. The la (3a) allowed range of 
Am 2 x = 8.2to;3(±o:s) x 10~ 5 eV 2 is decreased by a factor of 1.7 (5), but the allowed 
ranges of all other neutrino oscillation parameters and neutrino fluxes are not sig- 
nificantly changed. Maximal 9 12 mixing is disfavored at 5.8a and the bound on the 
mixing angle #13 is slightly improved to sin 2 6*13 < 0.048 at 3a. The predicted rate 
in a 7 Be neutrino-electron scattering experiment is 0.665 ± 0.015 (loo4o) °f the ra ^ e 
implied by the BP04 solar model in the absence of neutrino oscillations. The corre- 
sponding predictions for p— p and pep experiments are, respectively, 0.707^o;!h3(1o;q39) 
and 0.644^q'q^3(^q'q37). In order to clarify what measurements constrain which pa- 
rameters best, we also analyze the solar neutrino data separately and the reactor 
anti-neutrino data separately, both Before and After Neutrino 2004. We derive up- 
per limits to CPT violation in the weak sector by comparing reactor anti-neutrino 
oscillation parameters with neutrino oscillation parameters. We also show that the 
recent data disfavor at 91% CL a proposed non-standard interaction description of 
solar neutrino oscillations. We have verified that our results are insensitive (changes 
much less than la) to which of six approaches we use in analyzing the KamLAND 
data, which of the published 8 B neutrino energy spectra we adopt, and the precise 
value of the gallium solar neutrino event rate. 
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1. Introduction 

How have the recently released data from the the KamLAND reactor anti-neutrino 
experiment |lj and the revised average gallium solar neutrino rate |2|. |3|. f|, |[ improved 
our knowledge of neutrino properties and of solar neutrino fluxes? 

We concentrate in this paper on a Before-After comparison that is made pos- 
sible by the new data released at Neutrino 2004 (Paris, June 19-24, 2004). We 
determine how the new KamLAND and solar neutrino data affect our knowledge of 
the parameters that characterize solar neutrino oscillations [Am^, 6> 12 , #13, "^sterile] 
and the parameters that characterize solar energy generation and neutrino fluxes 
[£cno, <f>(p-p), 0( 7 Be), 0( 8 B), 0( 13 N), 0( 15 O), and 0( 17 F)]. 

In order to clarify which measurements constrain what quantities and by how 
much, we analyze the reactor data JIJ, ^|, [7j separately and the solar neutrino data || 
|3|, £|, ||, ||, |9], [10], [IT], [TJ, [L3|, [L4| separately We use six different approaches to analyzing 
the KamLAND data (see section |3.2| ) in order to assess the quantitative importance, 
or lack of importance, of different analysis procedures. 

The conventional wisdom is that a quantitative improvement in our knowledge of 
neutrino parameters and solar neutrino fluxes is all that we should expect. According 
to this view, the existing solar neutrino experiments have reached a level of maturity 
and precision at which new data from these operating experiments are expected to 
lead to refinements, but not revolutions. This conventional wisdom could be wrong 
and sub-dominant contributions due, for example, to non-standard interactions |T5, 



16fl , to sterile neutrinos [|T7| , |18| , or even to CPT violation |19[ could show up in the 
operating experiments. We investigate all these possibilities. 

We analyze the experimental data assuming that vacuum and matter neutrino 



oscillations |20], |2JJ occur among three active neutrino species (with the possibility 



also of oscillation to a sterile neutrino f22|). The techniques that we use in this 



analysis have been described previously in a series of papers, especially refs. [23, 



24, EH]. Many other groups have reported analyses of solar neutrino and reactor 



data, see ref . |26| , but our analysis is unique so far- we believe-in treating all of the 



solar neutrino fluxes as free parameters, subject only to the luminosity constraint, 
refs. [^7]. ^8j (i.e., essentially energy conservation). 

Our principal results are shown in figure [l| and figure ^] and in table ^] and table [3[ 
We begin by discussing in section B] the experimental data and the x 2 formula- 
tions we use for different applications. The data are summarized in table [I]. We then 
present in section |3] the results of our reactor-only analyses: the allowed regions in 
neutrino oscillation space that are compatible with the reactor data available Before 
Neutrino 2004 and the reactor data (notably the new KamLAND data [|l|) avail- 
able After Neutrino 2004. We also describe in section [^ the six different approaches 
we use in analyzing the KamLAND data and summarize the technical aspects of 
our KamLAND analyses. Next we present in section |] the results of our solar-only 
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analyses: the allowed regions in neutrino parameter space and neutrino fluxes that 
are compatible with all solar neutrino data available both Before and After Neu- 
trino 2004. Figure |1] summarizes the results of our Before- After reactor-only and 
solar-only comparisons. 

We present in section [| and in table |2| and table [3] the results of our global three 
neutrino analyses of solar neutrino experimental results and reactor anti-neutrino 
data. We give the best-estimates and the uncertainties for neutrino oscillation pa- 
rameters and for solar neutrino fluxes. We also determine in this section the upper 
bound on the sterile neutrino flux and on the luminosity of the Sun that is associated 
with the CNO nuclear fusion reactions. We compare in section [7] the allowed oscil- 
lation regions of rector anti-neutrinos with the allowed oscillation regions of solar 
neutrinos in order to establish an upper limit on CPT violation in the weak sector. 
We summarize and discuss our main results in section |[ In the Appendix, we present 
some details of the analysis involving #13. 



2. Experimental data and % 2 

We summarize in this section the experimental data we use and the x 2 distributions 
that we analyze. 

Table [J] summarizes the solar 0, 0, f|, H H 0] , reactor [I], || |7| , and atmo- 



spheric |3(J data used in our global analyses that are presented in section [5]. The 
number of data derived from each experiment are listed (in parentheses) in the sec- 
ond column of the table. In the third column, labelled Measured/SM, we list for each 
experiment the quantity Measured/SM, the measured total rate divided by the rate 
that is expected assuming the correctness of, as relevant, the standard solar model 
and the standard model of electroweak interactions (i.e., no neutrino oscillations or 
other non-standard physics). 

We calculate the global x 2 by fitting to all the available data, solar plus reactor. 
For the analysis of the upper bound on #13, we also include data from the K2K 
accelerator experiment and from atmospheric measurements, see eq. ( |A.5[ ). Formally, 



the global x can be written in the form p3L 34 



Xglobal - Xsolar( Am 21) #12, #13, {/b, fBe, fp-p, /cNo}) 

+ xLmLAND( A ™L 012, M • (2.1) 

Depending upon the case we consider, there can be as many as nine free param- 
eters in Xsoiar> including, Am^,^, 9 13 , / B , / Be , f P - P , and / C no (3 CNO fluxes, see 
below). The neutrino oscillation parameters Am^, 8 12 , #13 have their usual meaning. 
The reduced fluxes /b, /bg, fp-p-, an d /cno are defined as the true solar neutrino 
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Experiment 



Observable (# Data) 



Measured/SM 



Reference 



Chlorine 

SAGE+GALLEX/GNC-t 
Super-Kamiokande 
SNO (pure D 2 phase) 



SNO (salt phase) 



KamLAND 

CHOOZ 

K2K 

Atmospheric 



Average Rate (1) 
Average Rate (1) 
Zenith Spectrum (44) 
Day-night Spectrum (34) 



Average Rates (3) 



Spectrum (10) 

Spectrum (14) 

Spectrum (6) 

Zenith Angle Distributions (55) 



[CC]=0.30±0.03 

[CC]=0.52±0.03 

[ES]=0.406± 0.013 

[CC]=0.31±0.02 

[ES]=0.47±0.05 

[NC]=1.01 ± 0.13 

[CC]=0.28±0.02 

[ES]=0.38±0.05 

[NC]=0.90±0.08 

[CC]=0.69±0.06 

[CC] = 1.01 ±0.04 

[CC33(i/ M ) = 0.70±g;iS 

[0.5-1.0] 



t SAGE rate: 66.9 ± 3.9 ± 3.6 SNU §]; GALLEX/GNO rate: 69.3 ± 4.1 ± 3.6 SNU @, |. 



Table 1: Experimental data. We summarize the solar, reactor, accelerator, and 
atmospheric data used in our global analyses. Only experimental errors are included in 
the column labelled Result/SM. Here the notation SM corresponds to predictions of the 
Bahcall-Pinsonneault standard solar model (BP04) of ref. |3l]] and the standard model of 
electroweak interactions |32J (with no neutrino oscillations). The new average gallium rate 
is 68.1 ± 3.75 SNU (see ref. [Q]). The SNO rates (pure D2O phase) in the column labelled 
Result/SM are obtained from the published SNO spectral data by assuming that the shape 
of the 8 B neutrino spectrum is not affected by physics beyond the standard electroweak 
model. However, in our global analyses, we allow for spectral distortion. The SNO rates 
(salt phase) are not constrained to the 8 B shape |12| . The K2K and atmospheric data are 
used only in the analysis of #13, which is discussed in Appendix |A|. 



fluxes divided by the corresponding values of the fluxes predicted by the BP04 stan- 
dard solar model f3l| . We extend in section |5.2| the formalism to include sterile 
neutrinos. 

The function XxamLAND depends only on Am^, 6 I 12 and # 13 . 

We marginalize Xgiobai making use of the function Xchooz+atm+K2k(^i3) that was 



obtained following the analysis of ref. |35[ of atmospheric [f3(J , K2K accelerator |p£ 
and CHOOZ reactor |7|] data (see also, refs. |R| [37|]). We have not assumed, as is 



often done, a flat probability distribution for all values of 9^ below the CHOOZ 
bound. The fact that we take account of the actual experimental constraints on # 13 
decreases the estimated influence of 6*13 compared to what would have been obtained 
for a flat probability distribution. 
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3. Reactor data 

We compare in section [3J] the allowed oscillation regions for Am^, tan 2 # 12 , and 



sin 2 #13 that are determined from the first KamLAND results ||, together with the 
CHOOZ data JT], with the oscillation regions determined by including the recently 
released new KamLAND data [|TJ. This Before- After comparison is illustrated in 
the two right-hand panels of figure |l[ We also show in section |3.1| that the new 
KamLAND data more strongly disfavor a proposed [[15J non-standard description of 
solar neutrino oscillations. 



We describe in section |372| six different sets of assumptions that were used in 
analyzing the KamLAND data. We then discuss in section [3.3| some technical aspects 
of the analysis of the second release of KamLAND data. 

All of the neutrino properties and the solar neutrino fluxes that are determined 
in section [5] are robust with respect to the six different analysis approaches described 
in section 



3.1 Allowed regions: KamLAND reactor data 

In this subsection, we discuss briefly in section |3.1.1| the implications of the spec- 
tral distortion observed recently by the KamLAND collaboration, and then in sec- 



tion |3.1.2| present the best-fit values for neutrino parameters and their uncertainties, 
as well as the allowed contours obtained using the new KamLAND data. We show 
in section |3.1.3| that the new KamLAND data more strongly disfavor a previously 
proposed description of solar neutrino oscillations in terms of non-standard interac- 
tions. 

3.1.1 Spectral distortion 

The new KamLAND data |3J confirm the expected deficit of V e due to oscillations 
with parameters in the LMA region. More importantly, the new data show the 
expected distortion of the energy spectrum. In their new paper [[lj , the KamLAND 
collaboration report a goodness-of-fit test for a scaled no-oscillation energy spectrum 
with the normalization fitted to the data. They find a goodness-of-fit of only 0.1%. 
We confirm that the hypothesis of an undistorted scaled spectrum can fit the data 
with less than 0.2% probability. 

As a consequence, the 3<r region from the After KamLAND-only analysis shown 
in the lower right hand panel of figure |l| does not extend to mass values larger than 
Am^i = 2 x 10~ 4 eV 2 . For the now-excluded large Am^ values, the predicted 
spectral distortions are too small to fit the KamLAND data. 

3.1.2 KamLAND-only: best-fit values, uncertainties, and allowed regions 

Figure [I] compares the allowed regions for anti-neutrinos as determined by the Kam- 
LAND reactor experiment before Neutrino 2004 (upper right panel) with the allowed 
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Figure 1: Allowed oscillation parameters: Solar vs KamLAND. The two left 
panels show the 90%, 95%, 99%, and 3a allowed regions for oscillation parameters that are 
obtained by a global fit of all the available solar data [||, ||, |4|, 10, [t^, [14|, |5|, [t^]. The two 
right panels show the 90%, 95%, 99%, and 3<r allowed regions for oscillation parameters 
that are obtained by a global fit of all the reactor data from KamLAND and CHOOZ [0, [71]. 
The two upper (lower) panels correspond to the analysis of all data available before (after) 
the Neutrino 2004 conference, June 14-19, 2004 (Paris). The new KamLAND data Q are 
sufficiently precise that matter effects discernibly break the degeneracy between the two 
mirror vacuum solutions in the lower right panel. 



regions after Neutrino 2004 (lower right panel). The two panels in the right column 
of figure [i] represent the Before- After summary of the effect of the new KamLAND 
data released at Neutrino 2004. 

Before Neutrino 2004, the best-fit solutions for the reactor data were Am^ = 



7.1 x 10~ 5 eV 2 and tan 2 12 = 0.52 & 1.9 (see ref. pf). Within the statistical precision 
of the first KamLAND data, matter effects were too small to provide a meaningful 
discrimination between the two octants for 9\2. 

After Neutrino 2004, the best-fit values (x 2 = 12.5) is shown in the lower panel 
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of figure |T| and is: 

Am 2 21 = 8.3±g*^(±J-g) x lO^eV 2 , tan 2 12 = 0.36+°^(+^ 2 ) (reactor data : After). 

(3.1) 



As pointed out in Ref. [£3[], matter effects, although small, cannot be neglected in 
the analysis of precise KamLAND data |l| . In the analysis of the currently available 
data, matter effects break the degeneracy between the two octants of the mixing 
angle and induce an extra A% 2 = 0.2 for what would otherwise be (in the absence of 
matter effects) the mirror minimum at tan 2 # 12 = 2.7. 



We describe in section |3.2| six different approaches to analyzing the KamLAND 
data. The uncertainties shown in eq. (|3.1| ) are la (3a) errors for KamLAND analysis 
option number 3 of section jO| . 

The best-fit values and uncertainties of Am^ and tan 2 # 12 are essentially inde- 
pendent of the six analysis options for KamLAND data. The best-fit value for Am^ 
is the same to the numerical accuracy of eq. ( p.l|) for all six analysis options and 
the range of the la uncertainty varies by only about ±0.15a. The best-fit value of 
tan 2 #i2 varies by about ±0.2a and the range of the la uncertainty varies by ±0.1a. 

Our results are in good agreement with those obtained by the KamLAND collab- 
oration. They report a best fit point at Am^ = 8.3 x 10~ 5 eV 2 and tan 2 # 12 = 0.41, 
which is within the range of the best fit points obtained with our six analysis pro- 
cedures and almost identical to our preferred best-fit point. Comparing the results 
of our binned analysis with the results of the event-by-event maximum likelihood 
analysis of KamLAND, we find only two "barely visible" differences: (i) the lower 
"island" is allowed at a slightly lower CL in all of our binned analyses; (ii) the CL at 
which the "best-fit island" extends into maximal mixing is below 95% CL in Ref. [|1|], 
while in our binned analysis we find maximum mixing is slightly above or below 95% 
CL depending on the particular analysis option we adopt. 



The new KamLAND data 0, together with the CHOOZ 0, K2K |29|, and 
atmospheric p| results, lead to the following allowed range of sin 2 6*13, 



S m 2 6» 13 = o.oo5±8:8^(±8:8gi). (3.2) 



For the six different analysis options discussed in section |3T2| , the best-fit value of 
sin 2 # 13 varies by less than ±0.1a and the range of the la uncertainty varies by 
±0.1a. Note, however, that the best-fit value of sin 2 #13 is not significantly different 
from zero. 

3.1.3 Non-standard interaction 

Non-standard flavor- changing neutrino-matter interactions could potentially play a 
profound role in solar neutrino oscillations, even if the non-standard interactions are 



much weaker than standard weak interactions. In ref. [15|, Friedland, Lunardini, and 
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Pefia-Garay proposed a non-standard description of solar neutrino and reactor oscil- 
lations that expanded the allowed regions for neutrino oscillation parameters beyond 
what was allowed by standard interactions. The preferred oscillation parameters for 
this non-standard interaction are: 

^ m 2i = 1-5 x 10~ 5 eV 2 , tan 2 ^12 = 0.39 (non — standard interactions : ref. ||L5|| ). 

(3.3) 
After the new KamLAND measurements, this solution is disfavored at the 91% 
CL for 2 dof. This is a significant improvement over the first KamLAND results, 
which disfavored at 78% CL the non-standard solution of eq. (|3.3|). 

3.2 Six methods of analyzing KamLAND data 

We have analyzed the new KamLAND data with six different approaches. We first 
enumerate the six sets of assumptions that were used and then comment on the 
differences between the various assumptions. We provide additional technical details 



in the following subsection, section p75 . 

1. Poisson statistics, our normalization, 13 energy bins 

2. Poisson statistics, KamLAND normalization, 13 energy bins 

3. Poisson statistics, our normalization, 9+1 energy bins 

4. Poisson statistics, KamLAND normalization, 9+1 energy bins 

5. Gaussian statistics, our normalization, 9+1 energy bins 

6. Gaussian statistics, KamLAND normalization, 9+1 energy bins 

We prefer to combine the four highest energy bins of the KamLAND, which have 
only 6 events in total, in order to reduce the fluctuations and to make the analysis 
more stable. In order to verify that this additional binning does not affect the final 
results, we performed separate analyses with the full 13 published KamLAND energy 
bins (1 and 2 above) and with 9 + 1 energy bins (3, 4, 5, and 6). As we shall see in 
section |3.3| , our best-fit normalization for the number of observed events agrees with 
the best-fit normalization of KamLAND but the two normalizations differ slightly 
(by less than la). We have therefore performed analyses with our normalization 
(1, 3, and 5 above) and separately with the KamLAND normalization (2, 4, and 
6). Finally, we have compared, with the same energy binning and normalization, 
the results obtained with Poisson statistics (item 3) with the results obtained with 
Gaussian statistics (item 5). 

Of the six possibilities listed above, we prefer number 3. This option relies 
totally on our own calculations, so it is an independent check of the calculations 



of the KamLAND collaboration 0. Moreover, option 3 minimizes the effects of 
fluctuations due to low statistics bins. 

Fortunately, we shall show in section [5] that the globally-inferred results for 
neutrino parameters and solar neutrino fluxes are essentially independent of which 
one of the six options we choose. We have already seen in section |3J] that all six of 
the analysis options yield consistent results to an accuracy of much better than la. 

3.3 Some technical details: reactor anti-neutrino analysis 

The present analysis is based on data taken from 9 March, 2002 through 11 January, 
2004 |I[. We take account of corrections due to, among other things, the spallation 
cut and the detection efficiency of the tagged signal of electron antineutrinos (see 
ref. |TJ). We assume a time independent correction due to maintenance and bad 
runs and normalize our results after fiducial cuts to the KamLAND total exposure 
of 766.3 ton-year. 

We included in our calculations the time dependences due to the turn on/off of 
the different reactors in Japan. We have tracked the power of Japanese reactors on 
http : //www.fepc — atomic.jp/publicinfo/unten/index.html; the web page is owned 
by The Federation of Electric Power Companies of Japan. At present, this Web 
page tabulates the operational days of the 52 Japanese reactors up to June 2003. 
Furthermore, we have also been tracking the reactor power on a weekly basis since 
April 2003. This information allowed us to account for the time variations of the 
power-averaged antineutrino baseline, which can be as large as 20% (in good agree- 
ment with the results presented by KamLAND [l|]). Other time dependences like 
variations of the reactor composition could not be tracked, but have been shown to 
be small p8 |. We use the time averaged fuel compositions 235 U: 238 U: 239 Pu: 241 Pu 
= 0.568: 0.078 : 0.297: 0.057. Non- Japanese reactors contribute to the KamLAND 
signal less than 3% and their flux contribution is assumed time independent. 

With all this information, we find that, in the absence of oscillations, the expected 
number of antineutrino events above 2.6 MeV energy threshold is 381 which is in 
good agreement with the KamLAND estimate of 365.2 ± 23.7(syst). All of the solar 
neutrino parameters we infer from a global solution of the solar plus reactor data 
are, to high accuracy (much better than la) independent of which normalization we 



adopt (see items 1 and 2, items 3 and 4, and items 5 and 6 of section |3.2| and the 
discussion of the results all six analysis approaches in section ||]). 

We analyze the KamLAND energy spectrum by making a x 2 fit to their binned 
energy spectrum. The KamLAND spectrum contains a total of 13 energy bins above 
2.6 MeV, with only 6 events in the four highest energy bins. In order to reduce 
the fluctuations associated with the small number of events in these four bins, we 
combine for options 3, 4, 5, and 6 of section RT2] the data of these high energy events 



into a single bin with E > 6 MeV (containing 6 events). For this 10 bin analysis, we 
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compute the results assuming that the binned data is Poisson distributed, 



XKamLAND 



mm. 



10 



2(aR i 



th 



R 



cxp, 



+ 2K xp ln 



R, 



cxp 



aR th 



+ 



[a 



a, 



(3.4) 



HVS 



or that the binned data is Gaussian distributed, 



Xk. 



amLAND 



mm, 



10 



(aR\ h 



-^exp) 2 (« ~ 1) 



CT. 



stat,i 



cr. 



syst 



with cr 2 tatji = -R*- We a ^ so use a X 2 completely analogous to eq. 
all 13 energy bins (options 1 and 2 of section |3.2|) . 



(3.5) 



4[) when analyzing 



In eqs. (|3.4[ ) and (|3~5|), a is an absolute normalization constant and <7 syst = 6.5% 
is the total systematic uncertainty from several theoretical and experimental sources 
(see table I of ref. [[[]). In our binned analysis we have neglected the shape distortion 
errors. Using the presently available information we have been able to compute the 
shape distortion errors due to the uncertainties in the energy scale and reactor V e 
spectra and found them to be < 0.8% and < 0.5% respectively in any of the bins. 
We include in R\ h the expected number of events in the presence of oscillations, 
including the backgrounds from accidental coincidences (2.69 ± 0.02 events) and 
spallation sources (4.8 ± 0.9 events). The accidental background contributes to the 
event rate in the first bin while the spallation background is distributed among all 
the energy bins and peaks at E ~ 5.6 MeV. 

Our statistical analysis is different from the one of the KamLAND collaboration. 
First, we include the effect of #13 as described in Sec. |A|. Second, KamLAND col- 
laboration performs an unbinned maximum likelihood fit. Such an event-by-event 
likelihood analysis provides a more powerful tool to extract information from the data 
(see for instance Ref. ||39|). At the moment, only the KamLAND collaboration can 
perform an eventy-by-event maximum likelihood fit since to do so requires knowing 
the antineutrino energy (and time) for each event, which is not publicly available. 

In our previous studies |23|, |25], R3L EM , we used a calculational grid of 80 points 



per decade of Am^ and 80 points per decade of tan 2 # 12 . The previous grid is 
not sufficiently dense to take full account of the accuracy in the currently available 
neutrino data. Hence we are now using throughout the present paper a grid of 180 
points for each decade of Am^, 180 points for each decade of tan 2 6*12, and a step 
size of 0.00125 for sin 2 13 . 

As mentioned before, matter effects cannot be neglected in the present analysis 
of KamLAND data. They are most important in the lowest island and slightly favor 
the light versus the dark side of the mixing angle. To estimate the size of matter 
effects in the present analysis we define, F(matter vs vacuum), as the fractional 
difference in the event rate for the KamLAND detector calculated with and without 
including matter effects in the Earth. We find that, within the lcr(3cr) allowed 
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region of the analysis of the KamLAND energy spectrum, the maximum value of 
|_F(matter vs vacuum) | corresponds to 0.4% (2.3%). The maximum change in \ 2 due 
to including matter effects in the KamLAND-only analysis is 0.5 (1.4) at la (3a). 

4. Solar neutrino analysis 

In this section, we compare the allowed oscillation regions determined from all pre- 
vious solar neutrino experiments (chlorine, Kamiokande, SAGE, GALLEX/GNO, 
Super-Kamiokande, SNO) 0, | §, | g 0, @ with the oscillation regions 



determined by including the slightly revised average gallium rate released at Neu- 
trino 2004 [[| with the previously available data. 

We present in section |4.1| the main scientific results of this solar-only analysis . 



In section [4.2| , we describe some technical details of our analysis of the solar neutrino 
data. 

4.1 Allowed regions: solar neutrinos 

How much have the new solar neutrino data changed the allowed regions? The 
answer is 'imperceptibly', as the reader can easily see by comparing the upper left 
panel of figure [3] with the lower right panel of figure |]. We challenge even the most 
sharp eyed of our colleagues to discern the difference. 

Figure [I] shows the allowed regions for all solar neutrino experiments before 
Neutrino 2004 (upper left panel) with the allowed regions after Neutrino 2004 (lower 
left panel) . The two panels in the left column of figure [I] represent the Before- After 
summary of the effect of the new SNO data released at Neutrino 2004. 

The allowed regions for solar neutrino oscillations presented in figure [l] are some- 
what larger than the regions obtained by other authors f26j. The reason is that 
we have allowed all of the neutrino fluxes to be free parameters subject only to the 



luminosity constraint ||27|| , which is equivalent to energy conservation if light element 
fusion is the source of the solar luminosity. Most other groups |26[ incorporate in 
their analysis the solar neutrino fluxes and their uncertainties that are predicted by 



the standard solar model |31|, £|(J 



One can give good arguments for either including, or not including, the predicted 
solar model fluxes in the phenomenological analysis. The sound velocities measured 
from helioseismology are in excellent agreement with the standard solar model pre- 
dictions [IH | lj and the SNO measurement of the total 8 B neutrino flux [|l^, [l3| is 



also in agreement with the solar model predictions. These confirmations of the solar 
model justify the inclusion of the solar model predictions either as priors or as part 
of the x 2 analysis. 

We prefer instead to allow all of the solar neutrino fluxes to be free parameters in 
order to separate cleanly the astronomy from the neutrino physics. However, we have 
calculated the allowed ranges of neutrino oscillation parameters and neutrino fluxes 
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both ways, including the solar model predictions and letting all the neutrino fluxes 
be free parameters. Both methods yield similar-but not identical-results for Am^ 



and tan 2 ^12, although the method with free fluxes and the luminosity constraint 
yields a more accurate determination of the p — p solar neutrino flux |25| . 



4.2 Some technical details: solar neutrino analysis 



Details of our solar neutrino analyses have been described in previous papers [f23, 
P4] , j2"5"|| . The solar neutrino data we use are described in table ID. Solar data includes 
the Gallium (1 data point) and Chlorine (1 data point) radiochemical rates, the 
Super-Kamiokande zenith spectrum (44 bins), and SNO data previously reported 
for phase 1 and phase 2. The SNO data set available so far consists of the total 
day-night spectrum measured in the pure D 2 phase (34 data points), plus the 
total charged current (CC, 1 data point), electron scattering (ES, 1 data point), and 
neutral current (NC, 1 data point) rates measured in the salt phase |12], [U| [14]]. We 
use for the radiochemical experiments the neutrino absorption cross sections given 



in refs. g2], |3| 



We discuss in section |4.2.1| our treatment of solar neutrino fluxes. We discuss 
in section |4.2.2| how the choice of different available determinations of the shape of 
the 8 B neutrino energy spectrum affects the neutrino parameters and solar neutrino 
fluxes that are inferred using the existing solar neutrino and reactor anti-neutrino 
data. 

4.2.1 Treatment of solar neutrino fluxes 

Total neutrino fluxes are not required in our analysis and we only use the model 
fluxes to make dimensionless the neutrino flux output of our analysis. We express all 
neutrino fluxes determined by our phenomenological analysis of experimental data 
as ratios of the measured to the predicted (by the standard solar model BP04, |3l| ) 
neutrino fluxes. As a result of an obsessive sense for precision, we have used the most 
up-to-date electron and neutron densities and distributions of neutrino fluxes that 
are available on http: /www. sns.ias.edu/ jnb. We have checked, however, that none of 
our conclusions are affected significantly (<0.5%) by the particular choice of profiles 
we adopt. We have used solar models available at http:/www. sns.ias.edu/~jnb from 
2004, 2000, 1998, and 1995. Within the accuracy of the parameter determinations 
and the published grid sizes of the models, our inferences about neutrino parameters 
are unaffected by the choice of profile. All of the profiles are essentially unchanged 
to the accuracy of interest for solar neutrino work by the improvements in the mod- 
els over time (although the grid size has increased monotonically since 1962). We 
may regard the temperature and density profiles as experimentally confirmed be- 
cause of the remarkable agreement between standard solar model predictions and 
helioseismological measurements EO, ET . 
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All neutrino oscillation parameters and neutrino fluxes are treated as free pa- 
rameters. We obtain the allowed ranges of a particular parameter, solar or neutrino, 
by marginalizing the x 2 over an other parameters. In the x 2 analysis, the SNO and 
Super-Kamiokande sectors are correlated by the theoretical 8 B spectrum. 

4.2.2 Sensitivity to 8 B neutrino spectrum 

In our previous global studies of solar neutrino and reactor anti-neutrinos, we have 



used the Ortiz et al. 44] central values of the B neutrino energy spectrum and 



the Bahcall et al. errors |42[ on the energy spectrum. And, indeed, we used this 
prescription in all of the initial calculations in this paper. 

Recently, however, Winter et al. |EJ have redetermined the 8 B neutrino spectrum 
from new measurements of the a energy spectrum following the (5 + decay of 8 B. We 
have redone our three neutrino analyses replacing the Ortiz el al. central values of the 
8 B neutrino spectrum by the Winter et al. central values. We find that the Gallium 
and Chlorine rates are shifted downward by less than 0.1 a. For 8 B neutrinos, the 
charged current and neutral current rates on D 2 0, as well as the v — e scattering 
rate, are shifted downward by less than 0.1a, 0.3a, and 0.1a, respectively. The 
higher energy bins of the charged current and the electron scattering spectrum can 
be shifted downward by less than or equal to 5%. With the current state of the solar 
neutrino measurements, all of these changes are small compared to the uncertainties 
that result from a la shift in the experimental (electron) energy scale (up to 20% at 
la) and energy resolution errors (up to 16% at la). 

We have verified that the values of neutrino parameters and the values of neutrino 
fluxes derived in this paper from the different analyses are, to the accuracy we quote 
the numbers, independent of whether we use the Ortiz et al. 8 B neutrino energy 
spectrum or the Winter et al. energy spectrum. We have also verified that our 
results are unchanged to the quoted numerical accuracy if we use the Winter et 
al. error estimate on the shape of the energy spectrum or the more conservative 
Bahcall et al. error estimate. 

5. New global solution: solar plus reactor data 

In this section, we present the results of a global analysis of all the available solar and 
reactor data. We use the data summarized in table |l] and the total Xgiobai defined by 
eq. (pTl|). We marginalize Xgiobai with respect to 6> 13 by making use of the function 
Xchooz+atm+K2k(^i3) that is obtained from an analysis of atmospheric |3(J, K2K 
accelerator [E[J, and CHOOZ reactor M data (see Appendix for details). 

We begin in section |5.1| by presenting the solar neutrino oscillation parameters 
that are allowed Before and After Neutrino 2004 by the totality of existing neutrino 
data. We then summarize in section |5.3| the allowed ranges of solar neutrino fluxes. 



13 



Analysis Am| 1 (10 5 eV 2 ) tan 2 12 sin 2 13 sin 2 ?] 

Before 7.3±°; 4 (1™) 0.41+°;^ (+°; 22 ) < 0.019 (< 0.054) < 0.10 (< 0.38) 

After 8.2±g:l (lj;g) 0.39+°;^ (±°;g) < 0.015 (< 0.048) < 0.09 (< 0.35) 

Table 2: Global allowed regions. The table presents the best-fit oscillation param- 
eters and the la (3a) ranges determined from solar and reactor data available Before 
Neutrino 2004 (top row) and After Neutrino 2004 (bottom row). We have included the 
results of the marginalized x 2 distributions for #13 derived from the analysis of Atmo- 



spheric, K2K and CHOOZ data [36|. The quantity sin 2 77 characterizes the sterile fraction 
(see ref. [gjj, 25] and the discussion in section |5.2|) . The table presents the value obtained 
when all of the solar neutrino fluxes are treated as free parameters and the luminosity 
constraint is imposed; a slightly stronger limit is found if the solar model predictions of the 
neutrino fluxes are taken into account. The value of Xmin = 79.9(before) and 86.0(after) 
Neutrino 2004. 

5.1 Neutrino oscillation parameters: best-fit values, uncertainties, and 
independence of analysis method 

Table H gives the best-fit oscillation parameters and their uncertainties determined 
from a three-neutrino global fit to all the available solar and reactor data. We 
compare results obtained with the data available Before Neutrino 2004 (upper row) 
with the results obtained using also new data that became available only during or 
after Neutrino 2004 (lower row). 

For all of the neutrino oscillation parameters, the best-fit values obtained After 
Neutrino 2004 lie well within the quoted la uncertainties. 

The uncertainty for Amf^ represents the only dramatic improvement between 
the Before and the After analyses. The la (3a) uncertainty in Am^ is reduced by 
a factor of 1.7 (5.1). This reduction in the uncertainty of Am^ is almost entirely 
due to the observation of a statistically significant spectral energy distortion in the 
new KamLAND data IJ. For all the quantities shown in table ^| except for Am| 1? 
the After uncertainties have been reduced by < 25% with respect to the Before 
uncertainties. 

The neutrino parameters inferred from the global analysis are insensitive to which 
of the six analysis options discussed in section p.2| we use for the KamLAND data. 
The inferred values of Amfi and their uncertainties are identical for all six options 
to the numerical accuracy shown in table 0. The best-fit value for tan 2 9 12 (sin 2 # 13 ) 
varies by 0.3a (< 0.1a) for the six options and the la and 3a uncertainties vary by 
only 10% (< 10%). 

Figure |2] compares the allowed regions of solar neutrino oscillation parameters 
that were permitted before Neutrino 2004 with the allowed regions that are permitted 
after Neutrino 2004. The upper (~ 1.4 x 10 -4 eV 2 ) allowed region and the connecting 
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Figure 2: 'Before' and 'After' oscillation parameters: Solar plus reactor mea- 
surements. The figure shows the 90%, 95%, 99%, and 3<r allowed regions for oscillation 
parameters that are obtained by a global fit of the available solar || ||, [J, 1C , [0|, |l4|, || |1~2| 
and reactor data [jl], || Q. Left (Right) panel corresponds to analysis of data available 
before (after) Neutrino 2004 conference. The 'after' panels contain data from either ref. |l| 
(new KamLAND data) or ref. || (new gallium data). 



region between the upper and lower (~ 7 x 10~ 5 eV 2 ) allowed regions, all of which 
is permitted by the data existing before Neutrino 2004, is disfavored at more than 
3cr after the data presented at Neutrino 2004 are included. The elimination of the 
higher Am^ regions is due to the greater precision of the new KamLAND data (l[| . 
Maximal mixing is strongly excluded both before and after Neutrino 2004, 



tan 2 # 12 < 1.0 at 5.5ct(5.8cx) Before (After) 



(5.1) 



5.2 Sterile neutrinos 

We consider in this subsection the constraints on the admixture of sterile neutrinos. 
We concentrate on scenarios in which all Am 2 's but one are large enough to be 
averaged out in the solar and reactor neutrino experimental setups [22] 1 . The 
A-v models invoked to explain the LSND signal E6[ together with the solar and 



atmospheric neutrino data are an example of the kind of model we consider here. 

There are strong constraints of u e admixtures with large Am 2 's from the lack 
of observation of oscillations of reactor v e at short distances ]7|, [47], [48]]. These 
constraints imply that in the scenarios considered here solar and reactor oscillations 
can be understood as a single wavelength oscillation of v e into a state that is a linear 



x The effects of sterile neutrino admixtures with two different mass scales have been discussed 
for solar neutrinos in Ref. [[L7J . 
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Analysis pp 7 Be 8 B Lcno % 

Before Nu2004 1.02±° : °2 (±°;°g) 0.91±g;^ (±g;5?) 0.88±° : °t (±8£) 0.0+ 2 ; 8 (+ ™) 

After Nu2004 1.01±g;g| (gg) Lggg gffl) 0.87±g;gt (±g;gg) 0-0±g:g (gj) 

Table 3: Solar neutrino fluxes. The table presents the best-fit solar neutrino fluxes 



normalized to the BP04 solar model fluxes 1 31 ] . The la (3a) uncertainty ranges are deter- 



mined from all the solar and reactor data available Before (upper row) and After (lower 
row) the Neutrino 2004 conference. The last entry entry in the table is the fraction of 
the total solar luminosity that is associated with CNO nuclear fusion reactions. The fit to 
existing solar neutrino experiments prefers a CNO luminosity that is less than 0.1%, see 



ref. \2§. 



combination of active (v a ) and sterile (i/ s ) neutrino states, u e — > cosr]u a + smr)v a 
where rj is the parameter that describes the active-sterile admixture. The sterile 
contribution to the solar neutrino fluxes can also be parameterized in terms of sin 2 r\ 
or, alternatively, in terms of a derived parameter /b, sterile, the sterile fraction of the 
8 B neutrino flux (see discussion in ref |p3fl ). 



The 1(3)<T allowed range for the active-sterile admixture is 

sin 2 r? < 0.09(0.35) at l(3)a (5.2) 

for our global analysis of the existing solar plus reactor data. The fundamental pa- 
rameter describing the sterile fraction is sin 2 r\ (see ref |[23|). However, it is convenient 
to think in terms of a sterile fraction of the flux, /b, sterile , which is potentially ob- 
servable in the Super-Kamiokande and SNO experiments. This range corresponds 
to 

/ B ,steriie = 0.0^(±aoo) (solar + reactor) . (5.3) 

The limits on sterile neutrinos are not changed significantly by the additional 
data released at Neutrino 2004 (see the last column of Table |2|). 

5.3 Solar neutrino fluxes 

Table ^ presents the best-fit solar neutrino fluxes and their uncertainties that are 
obtained from a global fit to all the solar and reactor data. The fluxes are normalized 
to the predicted fluxes in the standard solar model |JT . 



The most remarkable fact about table y] is that the inferred ranges of the solar 
neutrino fluxes are very robust. The fluxes and their uncertainties are practically 
unchanged by the data released at Neutrino 2004. The p-p and 8 B solar neutrino 
fluxes are determined (la) to 2% and 5%, respectively, and the 7 Be and the CNO 
fluxes (see below) are very poorly determined. 

All of the estimates of the solar neutrino fluxes and their uncertainties are in- 
sensitive to which of the six KamLAND analysis options (see section 372) we employ. 
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o 9 +0.3 H-1.0\ 
°- z -0.3 K-0.8J 


o 9+O.3 /+1.0\ 
°' z -0.3 K-0.8J 


o 9+O.3 /+1.0^ 
°- z -0.3 K-0.8J 


U - OJ -0.04 l-0.ll/ 1 


n QO+0.05 /+0.19N 
U ' OJ -0.04 l-0.ll/ 1 


fl QO+0.05 /+0.21 
u -°°-0.05 1-0.11 


i n o+0.02 /+0.05\ 
i - uo -0.02 V-0.07J 


i n i+0.02 /+0.06\ 
i - ui -0.02 K-0.06J 


n qq+0.02 /+0.07 

u - JJ -0.02 1-0.06 


n O7+0.04/+0.09N 
U - 8 '-0.04Uo.llJ 


n 87 +0.04/+0.09\ 
U -0'-0.04Uo.llJ 


oo+0.04/+0.09n 
u - oo -0.04l-0.12/ 


n 9r+0.85/+1.37\ 
u - zo -0.25l-0.25J 


i n o+0.24 /+0.77\ 
i - uo -1.03 V-1.03J 


i 9Q+0.26 /+0.74 
-■-■ zy -0.57 1-1.29 



Gallium Rate (SNU) 63.3 ± 3.6 68.1 ±3.75 77.8 ± 5.0 

Am 2 21 (10- 5 eV 2 ) 

tan 2 6>i2 

p — p 
8 B 

7 R „ n OK+0.85/+1.37\ i no+0.24 /+0.77\ i 9Q +0.26 /+0.74n 

De u - zd -0. 251-0. 25J ± - uo -1.03 1-1.03J 1,zy -0.57 1-1.29^ 

Table 4: Analysis of solar and reactor data assuming different values of the event rate in 
gallium solar neutrino experiments. The table presents results for three different assumed 
gallium solar neutrino event rates; all the other solar and reactor data available prior to 
Neutrino 2004 were used unchanged in the analyses. The three cases considered are a low 
gallium rate (average of GNO and SAGE for the period 1998-2003), the average gallium 
rate (GALLEX/GNO and SAGE 1991-2003), and a high gallium rate (GALLEX and SAGE 
1991-1997). The solar neutrino fluxes are treated as free parameters, constrained by energy 
conservation (the luminosity constraint), and are normalized to the BP04 predicted fluxes. 

The best-fit p — p neutrino flux and the lcr and 3a p — p flux uncertainties are the 
same for all six options to the numerical accuracy of table |3|. The best-fit 8 B neu- 
trino flux varies by only 0.25a, the la uncertainty is the same in all cases, and the 3a 
range varies by 5% for all six options. The range of solutions found for the 7 Be solar 
neutrino flux using the six different analysis options is negligible compared to the 
uncertainties that arise from the experimental data and that are shown in table ^. 

The global analysis yields the following values for the CNO fluxes using data 
available before Neutrino 2004: 

13 N = .0+^ ; 1 5 O = 0> Q+2.4. 17p = Q+2.2 Before ^^ 

Using all the data available after Neutrino 2004, the CNO fluxes are: 

13 N = 0.0+^; 15 O = 0.0l; 2 ^ 17 F = 0t 2 1 After. (5.5) 

The global fit to all the data yields a best-fit CNO fractional contribution to 
the solar luminosity that is less than 0.1%, much less than the range, 0.8% to 1.6%, 
implied |3l| by standard solar model calculations. However, the constraint shown in 



table [| is consistent with the BP04 solar model prediction at la. 

5.4 How sensitive are the global results to the precise gallium event rate? 

The average rate of the two gallium experiments, GALLEX/GNO and SAGE, was 
slightly higher in the earlier period of data taking, 1991-1997, than it was in the 
later period, 1998-2003 |49], Q. In the earlier period, the average rate of the two 
experiments was (77.8 ± 5.0) SNU, while in the latter period, the average rate was 
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(63.3 ± 3.6) SNU gj|. The grand average for the period in question, 1991-2003, is 
(68.1 ±3.75) SNU. 

V. Gavrin |49j has raised the possibility that the difference between the observed 
rates in the earlier period and in the latter period could be due to a time-dependence 
of the lower-energy p — p and 7 Be fluxes. In this phenomenological interpretation, 
the hypothesized variation could not be observed in other solar neutrino experiments 
because of a lack of sensitivity to the p — p or 7 Be solar neutrinos or because of 
the long time scale and small amplitude of the assumed time dependence. Only 
further long-term measurements with low energy solar neutrino detectors can settle 
experimentally this question. We address here instead a related issue. 

We investigate in this subsection the extent to which the precise value of the 
gallium rate influences the inferred neutrino properties and solar neutrino fluxes 
that are obtained by a global solution to all the available solar and reactor data. 

Table |] presents the results of three different global analyses corresponding 
to the three different assumed gallium event rates: (77.8 ± 5.0) SNU (1991-1997), 
(68.1 ± 3.75) SNU (1991-2003), and (63.3 ± 3.6) SNU (1998-2003). All the other so- 
lar neutrino and reactor anti-neutrino data available prior to Neutrino 2004 were 



used unchanged in the analyses, which were performed as described in section [472 . 

The results obtained from the global analyses are remarkably robust with respect 
to the assumed gallium event rate. Within the range considered from 63 SNU to 78 
SNU, the inferred values for Am^ and tan 2 # 12 are essentially the same for all three 
choices of the gallium rate. The three best-constrained solar neutrino fluxes, p — p, 
8 B, and 7 Be, also change very little (changes less than ler)as the gallium rate is varied 
within the allowed range. 

6. Predicted rates for 7 Be, p — p, and pep experiments 

In this section, we summarize the predicted rates for future v — e scattering experi- 
ments with 7 Be, p — p, and pep solar neutrinos. 

Whether or not the 7 Be solar neutrino flux is treated as a free parameter af- 
fects strongly the uncertainty, and the best estimate, of the predicted rate in a 7 Be 
v — e scattering experiment. If we perform the global analysis of solar plus reac- 
tor experiments assuming that the BP004 calculations for the solar neutrino fluxes 
and their uncertainties are valid then the predicted rate in a 7 Be v — e scattering 
experiment |51] is, with lcr (3a) uncertainties: 



[ 7 Be] u _ e = 0.665 ±0.015 (l{j;{j^) [BP04 prediction]. (6.1) 

Here [ 7 Be]^_ e is the v — e scattering rate in units of the rate that would be expected 
if the BP04 solar model were exactly correct and neutrinos did not oscillate. Thus, 
if we use the solar model calculation of the 7 Be neutrino flux and its uncertainty, 
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then the predicted event rate for the 7 Be rate experiment has a precision of ±3%. If 
instead we treat all of the neutrino fluxes as free parameters, then the uncertainties 
are much larger: 

[ 7 Be]^_ e = 0.63+^l(+^) [all fluxes free + luminosity constraint] . (6.2) 



The result for v— e scattering that is given in eq. Q6.1|) and eq. (|6.2|) includes neutrinos 



from all the solar neutrino fluxes that produce recoil electrons with energies in the 
range 0.25-0.8 MeV. Most of the recoil electrons in the selected energy range are 
produced by 7 Be neutrinos, although there are small contributions from CNO, pep, 
and p — p neutrinos (whose fluxes we estimate by taking the values from the BP04 
solar model PI|). The best-estimates given in eq. ( |6.1| ) and eq. ( |6.2|) would be 



changed by much less than lcr, to 0.664 and 0.68, respectively, if we included only 
events caused by 7 Be neutrinos. 

The predictions given in eq. ( p-l\) and eq. ( |6.2j ) are robust. The results obtained 
here are well within the lcr quoted uncertainties of the pre-Neutrino 2004 predictions 
(see ref. [p5|). 

The corresponding predicted rates map — pu — e scattering experiment are, 
with la (3cr) uncertainties: 

[p-pU = 0.707+^3(^039) [BP04 prediction], (6.3) 

and 

[p - p ] v _ c = 0.716±g;^|(±g;g|2) [ a11 fl uxes free + luminosity constraint] . (6.4) 

The predicted rates in a pep v — e scattering experiment are: 

[pepL_ e = 0.644+Hl3(-ao3?) [BP04 prediction] , (6.5) 

and 

\pep] u _ e = 0.6h2tomi(toiti) [ a11 fluxes free + luminosity constraint] . (6.6) 

The predicted rates of the p — p and pep solar neutrino scattering experiments 
are robust |25| . The p — p and pep rates are affected by less than lcr by the choice 



of whether or not we use flux predictions from the BP04 solar model or allow the 
fluxes to vary as free parameters. 

Since the p — p, pep, and 7 Be solar neutrinos sample the survival probability at 
different energies (as reflected by the different values given in the predicted rates 
shown above), it is important to measure v — e scattering rates for all three fluxes. 
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7. CPT bound 

In this section, we compare the regions of allowed oscillation parameters for neutri- 
nos and anti-neutrinos in order to constrain the violation of CPT in the neutrino 



sector |52|, |53|, [54], |55], [19], [56]]. Our discussion follows the pre-KamLAND analysis of 
ref. [ID |, in which expected results from KamLAND were simulated. 

Figure [3] shows, in the left-hand panel, the 90%, 95%, 99%, and 3a allowed 
regions for oscillation parameters that are obtained by a global solution to all the 
available solar data and, in the right-hand panel, to all the available reactor data. 
The juxtaposition of purely solar and purely reactor data in the same figure allows 
a visual comparison of the constraints on neutrino and on anti-neutrino oscillation 
parameters. The allowed regions obtained by a global solution to all the available 
solar and reactor data are shown in figure ^. 



We characterize, following ref. [Tj|, the sensitivity of reactor and solar neutrino 



experiments to CPT violation by the quantity 

(ACPT) = 2 [g w "**"j (7.1) 

Here both K uu and Rpp are computed for the present KamLAND experimental set 
up, but using, respectively, values for (Ami, ®v) an d for (Ami, Op) within the allowed 
regions of the solar-only analysis (section |4.1| and the lower left-hand panel of figure [l]) 
and of the KamLAND-only analysis (section |3.1| and the lower right-hand panel 
of figure |l]) at a given CL. Then, (ACPT) is the number of events observed in 
KamLAND minus the number of events that are expected from the solar-only analysis 
if neutrinos and anti-neutrinos had exactly the same oscillation parameters, divided 
by the average. 

The maximum value of (ACPT) is 

(ACPT) < 0.52 (1.01) at la (3a). (7.2) 

Matter effects, which simulate CPT violation, contribute less than 0.01 (0.05) to 
(ACPT) for neutrino parameters in the la (3a) KamLAND-only allowed regions 
(lower right-hand panel of figure [I]) and 0.08 (0.09) for neutrino parameters in the 
Solar-only allowed regions (lower left-hand panel of figure [I]). 

This experimental upper limit on (ACPT) can be used to test arbitrary future 
conjectures of CPT violation. Following ref. []Tpj, we consider an effective interaction 



which has been discussed by Coleman and Glashow p2, |53|, and by Colladay and 



Kostelecky |54" |, that violates both Lorentz invariance and CPT invariance. The 

C(ACPT)) = vtbf 1[X vl (7.3) 



interaction is of the form 
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where a, (3 are flavor indices, L indicates that the neutrinos are left-handed, and b is 
a Hermitian matrix. We discuss the special case with rotational invariance in which 
bo and the mass-squared matrix are diagonalized by the same mixing angle. We also 
assume that there are only two interacting neutrinos (or anti-neutrinos) and follow 
previous authors in defining rj as the difference of the phases in the unitary matrices 
that diagonalize Am 2 and the CPT odd quantity 5b, which is the difference between 
the two eigenvalues of b . 

When the relative phase is rj = 0, the survival probabilities of neutrino and anti- 
neutrinos take on an especially simple form. In the case of reactor antineutrinos, 
the oscillations occur in vacuum to an excellent approximation (matter effects are 
negligible in the range of parameters allowed by KamLAND data, see details in 



ref. ||): 

'AmU 5bL 



P PP = cos 9 



13 



sin 2#i2 sin 



^21 J 



4E„ 



(7.4) 



In the case of solar neutrinos, sensitivity to 8b in the oscillation phase is lost because 
of the fast oscillations due to the AtttJi term. In both cases, solar and reactor 
neutrinos, the mixing in vacuum is unchanged. The standard survival probability 
and mixing in matter are valid with a modified ratio of matter to vacuum effects 
given by 

2y/2G F cos 2 9 13 n e E u 
P ~ Am 2 21 + 25bE v ' [ ' 

We have reanalyzed solar and reactor neutrino data with the modified probabil- 
ities that contain an extra free parameter 5b. We obtain the \ 2 distribution on 5b by 
marginalizing over all other variables. The resulting upper bound on the violation 
of CPT invariance is: 

5b < 0.6 (1.5) x l(T 20 GeV, rj = at lex (3a). (7.6) 



8. Summary and discussion 

We summarize and discuss in this section our main conclusions. 

The most important feature of the recent KamLAND results is the detection 
of the spectral energy distortion that was expected for oscillations with parameters 
within the Large Mixing Angle region. The distortion implies that the 3<r region 
from the KamLAND-only analysis no longer extends to masses larger than Am^ = 
2 x 1(T 4 eV 2 . 

Comparing the upper and lower right hand (KamLAND-only) panels of figure [I], 
one can see clearly the effect of the observed spectral distortion in eliminating pre- 
viously allowed regions with a large Am^. In the now-excluded large mass regions, 
the predicted spectral distortions are too small to fit the KamLAND data. 
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As a consequence of the observation of spectral distortion, the largest change in 
the globally allowed regions of neutrino parameters is the elimination of the larger 
Am^j regions in the global solution. This constriction of the globally allowed regions 
can be seen visually by comparing the left and right hand panels of figure |2|. The 
large Am^ regions were previously allowed in the global solution at 99% CL. 

The globally allowed range of AtttJi is reduced by a factor of 1.7 (5) at la(3a) 
compared to the previously allowed range (see figure |2] and table |2|) . The best-fit value 
of Amjj is increased by about 15% by the new data released at Neutrino 2004, which 
represents a 1.4a shift from the previous global minimum in the before Neutrino 2004 
analysis. 

The allowed ranges, as well as the best-fit values, of tan 2 9u and sin 2 #13 are, as 
is shown in table [| not affected significantly by the new data release. For example, 
the best-fit value of tan 2 # 12 is changed by only 0.4a and the la range of uncertainty 
is decreased by 10%. For sin 2 6*13, the 3a bound is decreased by 10%. 

The amount by which maximal mixing (tan #12 = l)is disfavored is 5.4a from 
just solar neutrino measurements (before and after Neutrino 2004) and 5.8a from a 
global solution of all the solar and reactor data available after Neutrino 2004. 

The upper limit to the parameter sin 2 77 that characterizes sterile neutrinos is 
very little affected by the recent data release, as is also shown table |[ Allowing all 
of the neutrino fluxes to be free parameters, the present upper bound on sin 2 r\ is 
equivalent to a maximum sterile fraction of 6% for the 8 B solar neutrino flux (see 
section |5~2j ). 

We conclude that the parameters which describe solar neutrino oscillations are 
robustly determined, if not yet as precise as we would like. This is the first major 
result of the present reanalysis. 

Our second major result is that the allowed ranges of the p — p and 8 B solar 
neutrino fluxes are robustly determined by the existing solar neutrino and reactor 
data. 

Table |3| shows that both the best-fit values and the la and 3a allowed ranges 
of all three of the major solar neutrino fluxes, the p — p, 7 Be, and 8 B neutrinos, are 
essentially unaffected by the data released at Neutrino 2004. The ±2 % precision 
with which the p — p neutrino flux is determined is due to the luminosity constraint 
as well as the neutrino measurements. 

The 7 Be and CNO neutrino fluxes are very poorly determined (both Before and 
After Neutrino 2004) (see table § and eq. ( |5.5|) ). The upper limit to the fraction of 
the solar luminosity that is associated with CNO fusion reactions remains about a 
factor of six above the value predicted by the BP04 solar model (see also table |3|). 

The rate of a 7 Be neutrino-electron scattering experiment is rather well predicted 
if the BP04 solar model calculations of the solar neutrino fluxes are assumed and is 
rather poorly predicted if all the solar neutrino fluxes are allowed to be free param- 
eters in fitting the existing solar and reactor data. If the solar model predictions are 
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adopted, then the predicted rate is 0.665 ± 0.015 (lo!o4o) °f the ra ^ e implied by the 
BP04 solar model in the absence of neutrino oscillations. The rate predicted if the 
neutrino fluxes are treated as free parameters is 0.63+o;24(-a57) of the solar model 
prediction without oscillations. 

The predicted rates for p — p and pep solar neutrino experiments are robust and 
essentially independent of whether one uses fluxes from the BP04 solar model or al- 
lows the fluxes to vary freely subject to the luminosity constraint. If we use the BP04 
fluxes, the predicted rate of a p — p experiment is [p — p]„_ c = 0.707io!oi3(-o!o39) 
and the predicted rate for a pep experiment is [pep] I/ _ c = 0.644+q;q;[3(^q;q37). 

Since the 7 Be, p — p, and pep fluxes sample the survival probability curve at dif- 
ferent energies (see the results and discussion in section ||), it is important to measure 
the v— e scattering of all three of these fluxes. According to equation (8) and figure 1 



of Ref. [25], all three of these fluxes should be in the region of neutrino oscillation 
space in which vacuum neutrino oscillations dominate over matter oscillations. 

CPT violation in the neutrino sector is limited by the comparison of the al- 
lowed neutrino oscillation regions (left hand panels of figure |l|) with the allowed 
anti- neutrino oscillation regions (right hand panels of figure [I]). This comparison 
leads to an accurate experimental limit. This limit is presented in both a model 
independent and a model dependent way in section |7]. 

We have verified that our results are insensitive (changes much less than la) to 
which of six approaches we use in analyzing the KamLAND data (section |3.2| ), which 
of the published 8 B neutrino energy spectra we adopt (section |4.2.2| j, and the precise 
value of the gallium solar neutrino event rate (section | 
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A. Analysis details regarding #13 

Table presents our limit on 13 . In this section, we describe some of the details 
of the method used to obtain this limit. We shall see that the #13 bound from 
the atmospheric and CHOOZ data implies that three neutrino effects are small for 
solar and reactor experiments. However these three neutrino effects are not totally 
negligible and they contribute to the final bound on #13. 

The minimum joint description of atmospheric, K2K, solar, and reactor data 
requires that all the three known neutrinos take part in the oscillation process. The 
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mixing parameters are encoded in the 3x3 lepton mixing matrix which can be 
conveniently parameterized in the standard form: 






U = c 23 S23 x 1 x -s 12 c 12 (A.l 



where Cy = cos% and Sij = sin^-. Note that the two Majorana phases are not 
included in the expression above since they do not affect neutrino oscillations. The 
angles Q^ can be taken without loss of generality to lie in the first quadrant, 0„ G 
[0,tt/2]. 

We know from the analysis of solar and atmospheric oscillations that the mass 
differences satisfy: 

Am 2 Q = A ™2i « |Arr&| = Am 2 tm . (A.2) 

Under this condition, the joint three-neutrino analysis simplifies. For the solar and 
KamLAND experiments, the oscillations with the atmospheric oscillation length are 
completely averaged out and the survival probability takes the form: 

Pi: = sin 4 9 l3 + cos 4 6 13 P% , (A.3) 

where in the Sun P^ is obtained with the modified density N e — > cos 2 9\ 3 N e . The 
analysis of solar data constrains three of the seven parameters: Am^, 12 and #13. 

For atmospheric and K2K neutrinos, the solar wavelength is too long to be 
relevant and the corresponding oscillating phase is negligible. As a consequence, 
the atmospheric data analysis restricts Am^ ~ Am^, #23, and #13. The mixing 
angle 13 is the only parameter common to both solar and atmospheric neutrino 
oscillations; it potentially allows for some mutual influence. The effect of #13 is to 
add a v^ — > v e contribution to the atmospheric oscillations. Finally, for the CHOOZ 
experiment, the solar wavelength is unobservable if Am^ < 8 x f0~ 4 eV 2 . The 
relevant survival probability oscillates with a wavelength determined by Am^ and 
an amplitude determined by $13. 

The above considerations imply that the global analysis factorizes as 

Xglobal = Xsolar+Kamland( Am 2D #12, #13) + Xatm+K2K+CHOOZ ( Am 31' ^23, #13) ■ (A.4) 

Thus the 3-u oscillation effects in the analysis of solar and KamLAND data are 
obtained from the study of: 

Xglobal I marg(Am 21 , 6>i 2 , #31) = Xsolar+Kamland(Ara 21 , #i 2 , #13)+Xatm+K2K+CHOOZ |marg(#13) • 

(A.5) 
We denote by Xatm+K2K+CHOozlmar g the x 2 analysis of atmospheric, K2K and CHOOZ 
data after marginalization over Am^ and Q 2 z- 
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The present strong bound on 6 13 is mostly determined by Xatm+K2K+CHOozUar</(#i3) 
and arises from the non-observation of V e oscillations at CHOOZ with an atmospheric 
wavelength. Here we have used the results from our updated three-neutrino oscil- 
lation analysis [[$J of atmospheric [57|, K2K [|2"9" |, and CHOOZ [0 data which we 



briefly summarize. 

The basic techniques employed can be found in ref. |3"5|, but we have updated the 



atmospheric neutrino analysis to include the results from the final Super-Kamiokande 
phase I analysis |57| . In particular we make use of the new three-dimensional fluxes 



from Honda |58[ as well as the improved interaction cross sections which agree better 
with the measurements performed with near detector in K2K. We have also improved 
our statistical analysis of the atmospheric data (see ref. |37| for details). 



As a result of the inclusion of all of these effects, the allowed region for Am^ 
is shifted to lower values and the 3a bound on #i 3 from Xatm+K2K+CHOOzlmar S (6 l i3) 
weakened slightly (from sin 2 6*i 3 < 0.056 to 0.059). Also the best-fit value of sin 2 # 13 
is not exactly at sin 2 # 13 = but rather at sin 2 # 13 = 0.006 (although sin 2 6 l 13 differs 
from 0.0 only by 0.5cr). This fact that sin 2 #i 3 is not exactly zero is due to the 
atmospheric neutrino data, in particular to the slight excess of sub-GeV e-like events 
which are better described with a non- vanishing value of 6 13 . 

The inclusion of the KamLAND and solar data further strengthens the #i 3 bound 
as is reflected in Table |[ Before the new KamLAND and solar results, the la (3a) 
bound on 6*i 3 from the global 3-v analysis was sin 2 6*i 3 < 0.019(0.054). This bound 
is now improved by about 10% to sin 2 # 13 < 0.016(0.048). 
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